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Abstract 

With large volumes of data available from LHC, it has become possible 
to study the multiplicity distributions for the various possible behaviours of 
the multiparticle production in collisions of relativistic heavy ion collisions, 
where a system of dense and hot partons has been created. In this context 
it is important and interesting as well to check how well the Monte Carlo 
generators can describe the properties or the behaviour of multiparticle pro¬ 
duction processes. One such possible behaviour is the self-similarity in the 
particle production, which can be studied with the intermittency studies and 
further with chaoticity/erraticity, in the heavy ion collisions. We analyse the 
behaviour of erraticity index in central Pb-Pb collisions at centre of mass 
energy of 2.76 TeV per nucleon using the AMPT monte carlo event gener¬ 
ator, following the recent proposal by R.C. Hwa and C.B. Yang, concerning 
the local multiplicity fluctuation study as a signature of critical hadronization 
in heavy-ion collisions. We report the values of erraticity index for the two 
versions of the model with default settings and their dependence on the size 
of the phase space region. Results presented here may serve as a reference 
sample for the experimental data from heavy ion collisions at these energies. 
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1 Introduction 


Dynamics of the initial processes, that is the distributions and the nature of inter¬ 
actions of quarks and gluons, in the heavy ion collisions affect the final distribution 
of the particles produced [1]. Of the various distributions, multiplicity distributions 
play fundamental role in extracting first hand information on the underlying particle 
production mechanism. If QGP is formed at these energies the QGP-hadron phase 
transition is expected to be accompanied by large local fluctuations in the number 
of produced particles in the regions of phase space [2- Thus the study of fluctua¬ 
tions in the multiplicity is an important tool to understand the dynamics of initial 
processes and consequently the processes of strong interactions, phase transition 
and also to understand correlations of QGP formation p]. 

A comprehensive theoretical model which can explain and give answers to all 
the complexities of the physics involved at high energy and densities, as is created 
in the heavy ion collisions, is still not available. A successful model focussed on one 
aspect of the problem may not say much about the other aspects, but at least should 
not contradict what is observed. The measures which are studied in the present 
work rely on the large bin multiplicities. At LHG energies multiplicities are high 
and it is possible to have detailed study of the local properties in ( 77 , (j)) space for 
narrow pr bins and thus to explore the dynamical properties of the system created 
in the heavy ion collisions. Thus as an initial attempt to understand the nature of 
global properties, as manifested in local fluctuations, here we develop and test the 
methodology and effectiveness of the analysis, analysing simulated events for Pb-Pb 
collisions at y/SN-j^= 2.76 TeV using A Multi-Phase Transport (AMPT) model. 

Study of charged particle multiplicity fluctuations is one of the sensitive probes 
to learn about the properties of the system produced in the heavy ion collisions. 
Factorial moments are one of the convenient tools for studying fluctuations in the 
particle production. The concept of factorial moments was first used by A. Bialas 
and R. Peschanski |4j to explain unexpectedly large local fluctuations in high mul¬ 
tiplicity events recorded by the JAGEE Gollaboratin. Advantage of studying fluctu¬ 
ations using factorial moments is that these filter out statistical fluctuations. The 
normalised factorial moment Fq is defined as 




(n!/(n- g)!) 
(n)'? 


( 1 ) 


where n is the number of particles in a bin of size in a d-dimensional space 
of observables and (...) is either vertical or horizontal averaging, q is the order of 
the moment and is a positive integer > 2. Then a power-law behaviour 


Fq{d) oc S-'f’” 


( 2 ) 


over a range of small 6 is referred to as intermittency. In terms of the number 
of bins M oc 1/(5, Eq.|^may be written as 


Fq{M) oc (3) 

where ipq is the intermittency index, a positive number. 
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Even if the scaling behaviour in Eq. is not satisfied, to a high degree of 
accuracy, Fq satisfies the power law behaviour 

Fq cx F^^ (4) 

This is referred to as F-scaling. In attempts to quantify systems undergoing 
second order phase transition, in Ginzburg-Landau (GL) theory [5], it is observed 
that 


/3, = (g-l)", = 1.304, (5) 

the scaling exponent, z/, is essentially idependent of the details of the GL pa¬ 
rameters. 

Factorial moments {Fq's) do not fully account for the fluctuations that the 
system exhibits. Vertically averaged horizontal moments, can gauge the spatial 
fluctuations, neglecting the event space fluctuations. On the other hand, horizon¬ 
tally averaged vertical moments lose information about spatial fluctuations and only 
measure the fluctuations from event-to-event. Erraticity analysis introduced in [6], 
where one finds moments of factorial moment distribution, takes into account the 
spatial as well as the event space fluctuations. It measures fluctuations of the spatial 
patterns and quantifies this in terms of an index named as erraticity index (/ig). In 
a recent work [7], is observed to be a measure sensitive to the dynamics of the 
particle production mechanism and hence to the different classes of quark-hadron 
phase transition. 

In erraticity analysis, event factorial moments are studied, defined for an 
event as 


fn{M) 

F^(M)= ® 

wherein. 

(6) 

fq{M) = {Umirim “ 1). {rim “ 9 + l))zi 

(7) 

where rim > 9 is the bin multiplicity of the bin, and (.. .)h 

over all bins such that for IVP cells 

is the average 

fq{M) = ^rnirim - 1)... (n^ - 9 -F 1) 

m—1 

(8) 

Now if Fq{M) fluctuates from event-to-event, then the deviation of Fq{M) from 
{Fq{M))y ((.. .)y is for averaging over all events) for each event can be quantified 
using order moments of the normalised order factorial (horizontal) moments 

that can be defined as 

Cp,,(M) = irqiM)). (9) 

where p is a positive real number and 


\F^{M)]P 

{Fq{M))P 

(10) 
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Figure 1: (a) Multiplicity distribution (b) pt (c) t] and (d) (j) distributions of 
charged particles generated in Pb-Pb collisions at ^/sJ^ =2.76 TeV using the 


DF and SM AMPT. 



Figure 2: dNch/drj versus Npart plot for DF and the SM AMPT, compared with 
the ALICE and CMS data 


To search for M-independent property of C'p q(M) one looks for a power-law be¬ 
haviour of Cp^q{M) in M, 

Cp,,(M) (11) 

this is referred to as erraticity |6j. \f t/jqip) is found to have a linear dependence on 
p, then erraticity index pq can be defined as 


_ dlpq{p) 

dp 


( 12 ) 


in the linear region so that it is independent of both M and p. fiq is a number 
that characterizes the fluctuations of spatial patterns from evevt-to-event. /j ,4 is ob¬ 
served |7] to be an effective measure to distinguish different criticality classes, viz., 
critical, quasicritical, pseudocritical and non-critical, having low value for critical 
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Pt 

window 

Default 

<N> 

String Melting 

<N> 

0.2 <pT < 0.3 

285.2 

434.8 

0.3 <PT< 0.4 

279.2 

355.5 

0.4 <PT< 0.5 

243.7 

271.6 

0.6 <pT < 0.7 

163.3 

155.5 

0.9 <pT < 1.0 

80.5 

66.1 


Table 1: Average Multiplicity of the Simulated Data sets analyzed in different 
Pt windows 



Figure 3: {X{rj), X{(p)) phase space of an event with M = 32 in DF and SM 
case 


hadronization compared to those having random hadronization. To a good approxi¬ 
mation, it is observed [7] that for the model with contraction owing to confinement, 
/r 4 (critical and quasicritical case) = 1.87±0.84 and for models without contraction 
/i 4 (pseudocritical and noncritical) = 4.65±0.06. These model values are suggestive 
of the significance of erraticity index to characterize dynamical processes. 


2 Data Analysed 

Charged particles in \r]\ < 0.8 and full azimuth, generated using two versions of 
A Multiphase Transport (AMPT) model [8] [9], in Pb-Pb collisions at -y/SNN = 
2.76 TeV are analysed. AMPT model is a hybrid model that includes both initial 
partonic and the final hadronic state interactions and transition between these two 
phases. This model addresses the non-equilibrium many body dynamics. Depending 
on the way the partons hadronize there are two versions, default (DF) and the string 
melting (SM). In the DF version partons are recombined with their parent strings 
when they stop interacting and the resulting strings are converted to hadrons using 
Lund String Fragmentation model. Whereas in the SM version, a quark coalescence 
model is used to obtain hadrons from the partons. 

We have generated 23424 DF and 19669 SM events with impact parameter, 
6 < 5, using the model parameters, a = 2.2, b = 0.5, p, = 1.8 and a = 0.47 . 
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Figure 4: Average bin multiplicity dependence on M for the five px bins, in 
case of DF and the SM AMPT model 


Multiplicity, px, pseudorapidity and (j) distributions of the simulated events is shown 
in the Figure[^ Charged particles generated in the \ri\ < 0.8 and full azimuth having 
px < 1.0 GeV/c in the small px bins of width 0.1 GeV/c are studied for the local 
multiplicity fluctuations in the spatial patterns. Five px bins are considered in the 
present analysis, as tabulated in the Table 1, along with the average multiplicity of 
the generated charged particles in the respective px bins. 

Though AMPT does not contain the dynamics of collective interactions that are 
responsible for critical behaviour, but it is a good model to test the effectiveness 
of the methodology of analysis for finding observable signal of quark-hadron phase 
transition (intermittency analysis) and the quantitative measure of critical behaviour 
of the system (erraticity analysis) at LHG energies. 


3 Analysis and Observations 


For an ‘e*^’ event, the order event factorial moment {F^{M)) as defined in Eq. 
([^ are determined so as to obtain a simple characterization of the spatial patterns 
in two dimensional (ry, (j)) space in narrow px windows. However we first obtain flat 
single particle density distribution using cumulative variable X{r]) and X{(f)) (lO) . 
which are defined as 


x{y) 


ly::: p^y^^y 


(13) 


here y is ry or cj), ymin and ymax denote respectively the minimum and maximum 
values of y interval considered, rj and (f) is mapped to X(?y) and X(())) between 0 
and 1 such that p{y) is the single particle rj or (p density. {X[r]), X{(p)) unit square 
of an event in a selected px window, is binned into a square matrix with bins 
where the maximum value that M can take depends on the multiplicity in the Xpx 
interval and the order parameter, so that the important part of the M dependence 
is captured. 
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Figure 5: P{F^) distributions for order moment q = 2 and g = 4 for DF and 

SM (0.3 < pt ^ 0.4). M values in multiples of 2 are shown only. 


To give a visualization of the binning in the (? 7 , </>) space in narrower bin, alego 
plot for an arbitrary event from DF AMPT data, in 0.2 < px < 0.3 window, with 
M = 32, is shown in Fig.|^ 

As value of M and px increases, the (ry, (p) space becomes empty, as is observed 
from Figure|^which shows the dependence of the average bin multiplicity ((n)) on 
M. Because of the denominator in Eq. ([^, a cluster of particles with multiplicity 
n > (7 in an event would produce a large value for F^{M) for that event. On the 
other hand, if the particles are evenly distributed, F^{M) would be smaller. Thus 
the spatial pattern of the event structure should be revealed in the distribution of 
F^{M) after collecting all events. 

Study of the event factorial moment distributions, {P{Fq)) reveal that the dis¬ 
tributions become wider as M increases and develop long tailsat higher q, especially 
at higher M values. Further the peaks of the distributions shift towards left with 
increase in M leading to decrease in (F®) (referred to as Fq hereafter) with M 
whereas the upper tails move towards right, at higher q values as M increases. It 
means that in small bins the average bin multiplicity (n) is so small that when there 
is a spike of particles in one such bin with n > q, the non-vanishing numerator in 
Eq. (1^ results in a large value for Fq{M) for that event. 

De^pendence of Fq on M can be studied in log-log plots as shown in Fig. for 
various px cuts. From the plots, it is observed that Fq's decrease as the bin size 
decreases or in other words, as M value increases. We observe for both the DF and 
SM in AMPT that relationship between Fq{M) and M is inverse of that in the Eq. 
([^; that is 

^AMPT(^) cx ^-<0 (14) 
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Figure 6: M Dependence of for DF as well as SM. 


Pt V- V- 


(Default) (String Melting) 


0.2 

< 

Pt 

< 

0.3 

1.738 

± 

0.008 

1.753 

± 

0.004 

0.3 

< 

Pt 

< 

0.4 

1.774 

± 

0.007 

1.793 

± 

0.005 

0.4 

< 

Pt 

< 

0.5 

1.758 

± 

0.006 

1.755 

± 

0.006 

0.6 

< 

Pt 

< 

0.7 

1.824 

± 

0.008 

1.869 

± 

0.016 

0.9 

< 

Pt 

< 

1.0 

1.778 

± 

0.013 

1.781 

± 

0.011 


Table 2: Scaling exponents for negative intermittency in the Default and String 
Melting versions of the AMPT Model 


Hence, with negative <p~ it is found that the charged particles generated by the 
default and the string melting version of the AMPT model exhibit negative inter¬ 
mittency. 

(141 suggests that (M) —?► 0 at large M and q, implying that in 

lere are too few rare high-multiplicity spikes anywhere in phase space. Eq. 


Eq 

AMPT t 


( [l^ is a quantification of the phenomenon exemplified by Fig. [^for one event, and 
is a mathematical characterization after averaging over many events. This same 
behaviour was observed in [7] for the events belonging to the non-critical class. 

We plot Fq versus F 2 in Fig. to check F-scaling. For each set of px bins 
linear fit has been performed to determine the value, jdq, the slope, as exemplified 
by the straight lines in Fig. |^(a). The dependence of jdq on {q — 1) is shown in 
Fig. which exhibits good linearity in the log-log plots. Thus we obtain a scaling 
exponent, denoted here as v_. In Tablej^are given the value of the negative scaling 
exponent for different windows and for both versions of the AMPT model studied 
here. Since the scaling that is there in Eq.j^is different to that we observe here for 
AMPT data, thus the scaling exponent obtained here cannot be compared with the 
V = 1.304 for the second order phase transition, as obtained from the GL theory. 


Since large fluctuations result in the high Fq tails of P{Fq), as exemplified 
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Fig. (b) and (d), it is advantageous to put more weight on the high Fq side in 
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Figure 7: F-Scaling for DF as well as SM mode of AMPT 




Figure 8: fiq versus [q — 1) plot for determination of the scaling exponents. 


averaging over P{Fq). That is just what the double moment Cp^q{M) does. We 
have determined C'p_q(M) for q = 2, 3, 4, 5 and p = 1.0, 1.25, 1.5, 1.75 and 2.0. 

The number of bins, M, takes on values from 2 to the maximum value possible 
while having reasonable {um) such that Fq ^ 0. To check whether Cp^q{M) follows 
the scaling behaviour with M, Cp^q is plotted against M. Fig. [^(a) to (d) shows 
respectively, for q — 2, Z, A and 5, the Cp^q versus M plot in the log-log scale 
for the window 0.6 < pt < 0.7 GeV/c, and for various values of p between 1 and 
2. As expected for all values of q, for p — 1.0, the Cp ^ = 1. For p > 1.0, Cp ^ 
increases with M and q values. Similar calculations are also done for the other px 
windows. In the high M region linear fits are performed for each q and p value so as 
to determine tpq{p). We see in Fig. 10 that for 0.6 < pt < 0.7 4’q{p) depends on 
p linearly for each q. Thus the erraticity indices defined in Eq. (12| are determined. 
Similar plots are obtained for the other px windows also and the values of pq are 
given in Table It can be seen from the table that as px value increases, the 
erraticity indices increase for both versions of the AMPT model. 

Comparing the pq values for the DF and SM data within the same window and 
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Figure 9: M dependence of Cp^q for the pr window 0.6 < pt < 0.7 GeV/c in 
case of DF and SM versions of the AMPT model 


for the same values of q, it is observed that pq has higher values for the DF version 
in comparison to SM for the pr windows below 0.6 GeV/c. That phenomenon 
is related to the average multiplicities of the two versions reversing their relative 
magnitudes at higher pT- However it is to be noted from Fig.jl^that the dependence 
of V’(j(p) on p is better distinguishable for the two versions of the AMPT for only q 
— 4. Coincidentally, as observed in [7], /i 4 seems to be a good measure to compare 
the erraticity indices of the different systems and data sets at these energies. 

We observe that the values of and /if^^ for all px windows are larger 
than those obtained for the critical data set in [7], on the same side as the non- 
critical case. We have found (p~ to be negative because P{Fq) broadens, as M 
increases, with {Fqc) shifting to the lower region of F®, thus resulting in negative 
intermittency. We did notice that the upper tails move to the right, suggesting 
the presence of some degree of clustering. To emphasize that part of P{Fq) we 
have taken higher p-power moments of (j)q{M), which suppress the lower side of 
Fq while boosting the upper side. The scaling properties of Cp^q{M) therefore 
deemphasize what leads to negative intermittency. Thus the erraticity indices pq 
reveal a different aspect of the fluctuation patterns than the scaling indices v~. 
Our study here has revealed interesting properties of scale-invariant fluctuations 
that should be compared to the real data. 


4 Summary 

Fluctuations in the spatial patterns of charged particles and their event-by-event 
fluctuations, as are present in the events generated using the default and string 
melting version of A MultiPhase Transport (AMPT) model are studied. This is a first 
attempt to study intermittency and erraticity at such high energies. It is observed 
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Figure 10: The ij^q dependence on the p for DF and SM AMPT in the 0.6 < 
Pt < 0.7 GeV/c window. 



AMPT 

(Default) 

AMPT 

(String Melting) 

Pt 



0.2 < Pt U.iJ 

0.3 <pT < 0.4 
0.4 <pT < 0.5 
0.6 <pT < 0.7 
0.9 <pt < 1-0 

0.043 ±0.002 

0.045 ± 0.002 
0.062 ±0.003 
0.154 ±0.008 
0.739 ±0.043 

0.016 ± 0.001 

0.027 ±0.001 
0.048 ± 0.002 
0.174 ±0.011 
1.014 ±0.064 




0.2 ^ Pt ^ 0.3 

0.3 <pt < 0.4 
0.4 <pT < 0.5 
0.6 <pT < 0.7 
0.9 < Pt 1.0 

0.901 ±0.081 

0.904 ±0.081 
1.304 ±0.118 
2.678 ±0.155 
4.502 ±0.147 

0.328 ± 0.025 

0.531 ± 0.043 
1.019 ±0.077 
2.832 ±0.137 
4.960 ±0.305 


P4 

Pa 

0.2 <pT< 0.3 
0.3 <pT < 0.4 
0.4 <pT < 0.5 
0.6 <pt < 0.7 
0.9 <pt < 1-0 

4.325 ±0.243 

4.532 ±0.234 
5.478 ±0.258 
5.640 ±0.203 
7.484 ±0.361 

2.481 ±0.183 

3.385 ±0.235 
3.935 ±0.021 
6.101 ±0.214 
7.359 ±0.305 


Ms 

Ms 

0.2 <pt< 0.3 

0.3 < Pt ^ 0.4 
0.4 <pT < 0.5 
0.6 <pT < 0.7 
0.9 < Pt ^ 1.0 

6.202 ± 0.302 

6.150 ±0.175 
7.396 ±0.437 
9.107 ±0.693 
8.643 ±0.537 

5.143 ± 0.022 

5.745 ±0.312 
6.159 ±0.280 
8.360 ±0.533 
7.655 ±0.358 


11 














that as the bin size decreases, the factorial moments decrease. This behaviour is in 
contrast to usual properties of intermittency observed at lower energies indicating 
that events with localization of even moderate multiplicities in the small bins, at low 
Px are not present in the AMPT. Further the erraticity analysis of the model shows 
that the systems generated in it is not near criticality. The pg values determined 
here give the quantification of the event-by-event fluctiations in the spatial patterns 
of the charged particles in the midrapidity region, which can be used effectively to 
compare with other models. More importantly comparison of the values with that 
from the LHC would help to get a better understanding of the particle production 
mechanism at high energies. 
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